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Abstract 

Canonical formulation of quantum field theory on the Light Front (LF) is re- 
viewed. The problem of constructing the LF Hamiltonian which gives the the- 
ory equivalent to original Lorentz and gauge invariant one is considered. We de- 
scribe possible ways of solving this problem: (a) the limiting transition from the 
equal-time Hamiltonian in a fastly moving Lorentz frame to LF Hamiltonian, (b) 
the direct comparison of LF perturbation theory in coupling constant and usual 
Lorentz-covariant Feynman perturbation theory. Gauge invariant regularization of 
LF Hamiltonian via introducing a lattice in transverse coordinates and imposing 
periodic boundary conditions in LF coordinate x~ for gauge fields on the interval 
\x~\ < L is considered. We find that LF canonical formalism for this regularization 
avoid usual most complicated constraints connecting zero and nonzero modes of 
gauge fields. 



*Published in the book: "Quantum Theory in honour of Vladimir A. Fock", Part 1. Unesco, St. Pe- 
tersburg University, Euro-Asian Physical Society, 1998. P. 38-97. 
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1 Introduction 



V. A. Fock has elaborated a beautiful method to describe state vectors in quantum field 
theory. The corresponding vector space is now called the Fock space This method 
plays very important role in quantum field theory in Light-Front (LF) coordinates [0: 
x ± = (l/\f2)(x° ± x 3 ), x 1 , x 2 , where Lorentz coordinates. The x + plays 

the role of time, and canonical quantization is carried out on a hypersurface x + = const. 
The advantage of this scheme is connected with the positivity of the momentum P 
(translation operator along x~ axis), which becomes quadratic in fields on the LF. As a 
consequence the lowest eigenstate of the operator P_ is both physical vacuum and the 
"mathematical" vacuum of perturbation theory ||. Using Fock space over this vacuum 
one can solve stationary Schroedinger equation with Hamiltonian P + (translation operator 
along x + axis) to find the spectrum of bound states. The problem of describing the 
physical vacuum, very complicated in usual formulation with Lorentz coordinates, does 
not appear here. Such approach, based on solving Schrodinger equation on the LF, is 
called LF Hamiltonian approach. It attracts attention for a long time as a possible mean 
for solving Quantum Field Theory problems. 

While giving essential advantages, the application of LF coordinates in Quantum Field 
Theory leads to some difficulties. The hyperplane x + = const is a characteristic surface for 
relativistic differential field equations. It is not evident without additional investigation 
that quantization on such hypersurface generates a theory equivalent to one quantized in 
the usual way in Lorentz coordinates ||, ||, |6|, [7|, || || [10], [II], |I|, [13| . This is in particular es- 
sential because of the special divergences at p_ = appearing in LF quantization scheme. 
Beside of usual ultraviolet regularization one has to apply special regularization of such 
divergences. We will consider the following simplest prescription of such regularization: 

(a) cutoff of momenta p_ 

|p_|>e, s>0; (1.1) 

(b) cutoff of the x~ 

-L<x~<L. (1.2) 

with periodic boundary conditions in x~ for all fields. 

The regularization (b) discretizes the spectrum of the operator P (p_ = irn/L, where 
n is an integer). This formulation is called sometime "Discretized LF Quantization" HHf . 
Fourier components of fields, corresponding to p_ = (and usually called "zero modes") 
turn out to be dependent variables and must be expressed in terms of nonzero modes via 
solving constraint equations (constraints). These constraints are usually very complicated, 
and solving of them is a difficult problem. 

The prescriptions of regularization of divergences at p- = described above are con- 
venient for Hamiltonian approach, but both of them break Lorentz invariance and the 
prescription (a) breaks also the gauge invariance. Therefore the equivalence of LF and 
original Lorentz (and gauge) invariant formulation can be broken even in the limit of 
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removed cutoff. To avoid this inequivalence some modification of usual renormalization 
procedure may be necessary, see for exsample |L5| and [0 . 



The problem of constructing the LF Hamiltonian which gives a theory equivalent to 
original Lorentz and gauge invariant one turned out to be rather difficult. Now it is solved 
only for nongauge field theories We will describe possible approaches to this problem. 

In Sect. 2 we give basic relations of quantum field theory in LF coordinates. In 
Sect. 3 we consider the limiting transition from fast moving Lorentz frame to the LF. 
This transition relates the formulations in Lorentz and LF coordinates [T2|]. In Sect. 4 we 
investigate the relation between LF perturbation theory in coupling constant and usual 
Lorentz-covariant Feynman perturbation theory. For Yukawa model this investigation 
shows how to construct LF Hamiltonian giving the theory perturbatively equivalent to 
original one For gauge theories the methods developed in Sect. 3 and Sect. 4 do 



not give the required LF Hamiltonian due to specific difficulties. In Sect. 5 we consider 
gauge invariant ultraviolet regularization of LF Hamiltonian via introducing a lattice in 

x 2 and taking complex matrix variables for gauge fields on 



transverse coordinates x 1 
links of the lattice 
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We find that LF canonical formalism for gauge theories with 
this regularization avoid usual most complicated constraints connecting zero and nonzero 
modes. 



2 Formal canonical quantization of Field Theory 

on the Light Front and the problem of bound states. 

In order to find the bound state spectrum in some field theory quantized on the LF the 
following system of equations is usually solved: 

p + l*> = p;i*>, (2.i) 

P_\V) =P'_\V), (2.2) 

Pj_|tt) = 0, (2.3) 
where P± = {Pi, P 2 }. The mass of bound state is equal to 

m = ^2P' + P'_. (2.4) 
It was taken into account that nonzero components of metric tensor in LF coordinates 

are 

9-+ = g+- = 1, 9n = 922 = -1. (2.5) 

The operators P_, P± are quadratic in fields, and the solution of equations (|2.2fl, fl2.3|) is 
not difficult. The problem is in solving the Schroedinger equation (|2.1| ). Physical vacuum 
\Q) is lowest eigenstate of the operator P_, and 

P + |fi) = 0, (2.6) 
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P_|0) = 0, 



(2.7) 



P±\Q) = 0. (21 



To fulfil equations (|2.6|), ( 2T7|) one should subtract , if it is necessary, corresponding 



renormalizing constants from the operators P+, P_. The plays simultaneously the 
role of mathematical vacuum of Fock space. A solution of Schrodinger eqn. Q2.1Q 
belong to this space. 

Expressions for the operators P + , P_ can be obtained by canonical quantization on 
the LF. Let us describe the procedure of such quantization via some examples, without 
analyzing so far the question about the equivalence of appearing theory and the orig- 
inal Lorentz-covariant one. It is assumed that in addition to explicit regularization of 
divergences at p_ = necessary ultraviolet regularization is implied. 

2.1 Scalar selfinteracting field in (l+l)-dimensional 
space-time. 

Peculiarities of LF quantization are well seen even in this simple example. We have only 
LF coordinates x + , x~ . The Lagrangian is equal to 

L = J dx~ (d^ip - -m 2 <p 2 - V) , (2-9) 



or 



L = J dx~ (^d+ipd^ip - ^mV - V 4 ) , (2-10) 

The "time" derivative d+ip enters into this Lagrangian only linearly. For the transition 
to canonical theory it is sufficient to rewrite the expression / dx~ \d+ipd-ip in standard 
form. To achieve this let us take the Fourier decomposition 

oo 

(f(x~) = {2ix)~^ J dk |2fc|~3 (a(k) exp(—ikx~) + a + (k) exp(ikx~)^ , (2.11) 



where k = k-, ip(x ) = f{x + , x ), a(k) = a(x + , k). The Lagrangian ( 2.10 ) takes the form 



L = I dk (* + (k)a(k)-a(k)a+(k)\ _ ^ ^ 



o 



where a = da/dx + and 

H = J dx~ ^mV 2 + Ay? 4 ^ . (2.13) 
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Here we have used the equality 



dk I dk'5{k + k')k' {a(k)a{k') - a + {k)d + {k')) = 0. (2.14) 



o o 



It is implied that the function <£>(x) in fl2.13|) is expressed in terms of a + (k) and a(k) 



with the help of formulae ( [2.10|) . "Time" derivatives d(k), d + (k) enter into Lagrangian 
L in a form standard for canonical theory. Therefore one can interpret after quantiza- 
tion the a + (k) and a(k) as creation and annihilation operators satisfying the following 
commutation relations at fixed x + and k > 0, k' > 0: 

[a(k), a + (k')] = 5{k - k'), [a(k), a(k')} = 0, (2.15) 

It is also seen that the H is LF Hamiltonian, i.e. H = P + . 
We have also the formulae 



P M = y dx-T^, (2.16) 

where the energy-momentum tensor T UfM is equal to 

Tvp, = dycpd^ip - g^C. (2.17) 

Via this relation one can reproduce the expression (|2.13|) for P + = H , and obtain the 
equality 



P_ 



oo 

dx~ {d^f = -Jdkk (a + {k)a{k) + a{k)a + {k)) . (2.18) 
o 



The lowest eigen state of the operator P_ is the physical vacuum \Q) for which 

a(k)\tt) = (2.19) 

at any k. It is seen that vacuum expectation values (Q\P-\Q), (Q\P + \Q) are infinite. The 
renormalization can be got by taking normal ordered forms : P + :, : P_ : with respect to 
the operators a + , a (the symbol :: means as usual that operators a + stand everywhere 
before of operators a). Normal ordering of the Ay? 4 term in the Hamiltonian P + leads also 
to renormalization of the mass. In the following, writing P + , P_, we mean the expressions 
: P + :, : P_ :, satisfying conditions ( |2.6| ), (|2~7j). Normal ordering of the operators P + , P 
allows to avoid all ultraviolet divergences in this simple model. 

2.2 Theory of interacting scalar and fermion fields in 
(3+1)- dimensional space-time (Yukawa model). 

The Lagrangian of the model is 

L = J d 2 x x dx- (V (i'fdn — M) ip + )^d^d^ - -^m V - 9^<P ~ ^V 3 - V 4 ) , (2-20) 
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where M is the fermion mass, m is the boson mass, ip = f = V 9 " 1 "! 9i \ A' are 

coupling constants. Here and so on fi, v, . . . — +, — , 1,2; i,k, ... — 1, 2; or 1 = (x l , x 2 ) . For 
Dirac's 7-matrices we use: 

( -%I \ 3 / il \ ± ( -to 1 - \ , 001 . 

where J is a unit 2x2 matrix, <j l = {a 1 , a 2 }, a 1 are Pauli matrices. 
We introduce 2-component spinors \, £, writing 



X 
t 



The Lagrangian L can be written in the form 

L = J d 2 x ± dx~ (iV2 X + d+x + iV2^ + d^ + (a% -Af)x + H - c -) + 
+d+(pd-<p - \d iV d iV - \m\ 2 - ig^ U + X ~ X + d ~ AV 3 - A^ 4 ) , (2.22) 



where H.c. means Hermitian conjugation. The variation of this Lagrangian with respect 
to x + leads to the equation 

V2d^ = - (<r% -M)x + 9fX (2.23) 



This equation does not contain derivatives in x + and therefore is a constraint. One 
should solve it with respect to £ and substitute the result into the Lagrangian. In doing 
this we must invert the operator d- which becomes an operator of multiplication ik- after 
Fourier transformation. Inverse operator (ikJ)^ 1 has singularity at k^ = 0. To avoid this 
singularity we introduce the regularization ( |1.1|) . For any function f(x~) = f(x + ,x~,x ± ) 
we define Fourier transform 

f(x~) = ^=Jdk. e ik -*-f(k.), (2.24) 



where /(&_) = f(x + ,k_,x ± ), and put 



[fix')) = — L [ dk.e ik ~ x ~ /(*_), |jfe_| > e > 0. (2.25) 



/2tt 

We insert into the Lagrangian ( |2.22| ) the variables [x], [£], [f] instead of x, 
X + 1 £5 £ + > V 9 an d obtain the constraint equation 



V2d4£] = - (a% - M) \x\+g[[<pM\ (2.26) 
instead of fl2.23| ). Its solution is 

[£] = ±fl\-(<r%-M) [x}+9[[M), (2-27) 
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where the operator d_ l is completely defined by the condition 

dZ 1 [f]= \dZ l [f] . 



(2.28) 



After Fourier transformation the operator dz 1 is replaced by (ikJ)~ . 

Substituting the expression ( p. 27] ) into the Lagrangian (where all fields x, X + , • • • are 
replaced with [x], [x + ]> . . .) we come to the result: 



L = J d 2 x ± dx~ UV2 [x + \ d + [x] + <9_[y?]<9 + [^] + 
+-^( (a% - M) [x] - g[[ V ][x}]) + {-id-)' 1 ( (a k d k - M) [ X ] - *][%]])- 

-\diWM - ^™ 2 M 2 - Avf - AM 4 ), 



(2.29) 



As in Sect. 2a time derivatives <9 + [x]> d + [ip} enter into the Lagrangian ( |2.29| ) linearly. 
Therefore to go to canonical formalism it is sufficient to find a standard form for the 
expression 

iV2[x + )d + [x]+d-[<P + ]dM 
(before quantization the quantities x + 1 X are elements of Grassman algebra). We write 



[np(x~)] = {2tt)- 1/2 I dk_ (2k_)- 1/2 {a{k_) exp{-ik_x~) + a + {k_) exp{ik_x')), (2.30) 

OO 

lXr(x-)] = (2n)- 1/2 2- 1/4 J dk-{b r {k~) exp(-ik-X~) + c£(A;_) exp(iJfe_aT)), (2.31) 

e 

oo 

[x+(aT)] = (2tt)- 1 / 2 2- 1 / 4 f dk„{c r {k„) exp{-ik_x~) + &+(&_) exp(ik-x~)), (2.32) 



where 

= [(f(x + , x~, x" 1 )], a{kJ) = a(x + , k-,x x ) 
et cetera, r = 1,2. The Lagrangian Q2.29Q takes the form 

oo 

d 2 x ± J dk_ (a(k„)d + (k_) - a + (k_)d(k_)- 

£ 

-b r (k-)b+(k.) - &+(fc_)6 r (fc_) - c r (A;_)c+(£;_) - c+{k^)c r {k_)) - H, 

where 



(2.33) 



- ? =( - Mj bd - #[[</?] [x]]J HS-)-^^ - mj m - *][x]]J + 



(2.34) 
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It is assumed that the quantities [x], [f] m the formulae ( |2.34| ) are expressed in 

terms of b, b+, c, c+ a, a + with the help of flE50D , ( fE3ID , ( ^32|) . 



It follows from (|2.33|) that a + , a, 6 + , 6, c + , c play a role of creation and annihilation 
operators. After quantization they satisfy the commutation relations (at x + =const): 

[a(k_, x x ), a + (k'_, = 5{k_ - k'_)8 2 {x L - x f± ), (2.35) 

[&(*;_, ar 1 ), &+(#_, x ,± )]+ = <5(A;_ - ^)5 2 (x ± - x f± ), (2.36) 

[c(Jfe_, ar 1 ), c + (fcl, x 7± )]+ = <S(A;_ - jfe / _)«J 2 (a:- L - x ,x ), (2.37) 

where [a;, y]± = xy ± yx. The remaining (anti) commutators are equal to zero. It is seen 
that the quantity H is LF Hamiltonian (H = P+). The operator of the momentum P is 
equal to 

1 °° 

P-= J d 2 xdx~T__ = -J (fx 1 - J dk_ k_(a + (k^)a(k_) + a(k_)a + (k_) + 

e 

+6+(fc_)6 r (A;_) - b r (k^)b+(k_) + c+{kJ)c r {kJ) - Cr{k-)c+{k-)). (2.38) 

The quantities P + = H and P_ should be normally ordered with respect to creation 
and annihilation operators. The lowest eigenstate of the momentum P_ is the physical 
vacuum. It is defined by conditions 

a(fc_,x ± )|fi) = 0, 6(^,^)1^) = 0, cik^x^lQ) = 0, Vi 1 ,^. (2.39) 



The equalities fl2.6|) , ( |2.7|) become true after normal ordering of the operators P + and P_. 
For the P_ it is seen from the formulae ( 2.38|) , and for the P + it follows from the fact 



that every term of P + contains a 5-function of difference between the sum of momenta k- 
of creation operators and the sum of momenta fc_ of annihilation operators. Due to the 
positivity of all momenta in our regularization scheme every term of the P + contains 
at least one annihilation operator. Therefore for normally ordered P + we have P+|0) = 0. 

The model under consideration requires ultraviolet regularization. It can be done in 
different ways. We consider this question together with the renormalization problem in 
Sect. 3 and 4. 

2.3 The U(N)-theory of pure gauge fields. 

We consider the U(N) rather than the SU(N) theory because it is technically more simple. 
The transition to the SU(N) can be done easily. Gauge field is described by Hermitian 
matrices 

A^x) = A+(x) = {A*{x)} , (2.40) 



8 



where \i = +,—,1,2; i,j = 1,2, . . . , N. Let us assume, that for the indexes i,j and 
analogous the usual rule of summation on repeated indexes is not used, and where it is 
necessary the sign of a sum is indicated. Field strengths tensor is 



d^Ay - d v A^ - ig[A^, A L 



(2.41) 



and gauge transformation has the form 



(2.42) 



An -> A\ L = U+A^U + -U+d^U, U + U = I. 

To escape a breakdown of gauge invariance we apply the regularization of the type (|L~2] 
with periodic boundary conditions 



Afji(x , L, x 



A^x + ,L,x x ), 



(2.43) 



on the interval — L < x~ < L. All Fourier modes of A^(x) in x~ must be kept, including 
zero modes (at fc_ =0). 

The Lagrangian has the form 



L 

L = — - J (fx 1 - J dx-Ti (F^F^) 



(2.44) 



or 



d 2 x A 



<lv \ r I F + _ + 2F_ k F +k — -Fkk'Fkk' ) ■ 



(2.45) 



where k, k' = 1, 2. Time derivatives d + Ak are present only in the term 

Tr (2(d^A k - d k A^ - ig[A_, A + ])d + A k ) . 

This expression can be put in standard canonical form only after fixing the gauge in a 
special form of the type A^ = because then the term above becames similar to that 
of scalar field theory Tr (2(d-A k )d + Ak) . However not every field, periodic in x~ , can be 
transformed to the A_ = gauge. Indeed, the loop integral 



F(x + ,x ± ) = Tr < 



P exp 



i 



dx A_(x + ,x ,x J 



(2.46) 



where the symbol 'P' means the ordering of operators along the x~ , is a gauge invariant 
quantity. If this integral is not equal to iV for some field then the gauge A_ = is not 
possible for this field. Therefore we choose more weak gauge condition 



d_A l 



0. 



(2.47) 
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and put 



A il L{x) =v i (x + ,x ± ). 



(2.48) 



The gauge (|2.47|) breaks not only the local gauge invariance but unlike the A_ = gauge, 
also global gauge invariance (it remains only the abelian subgroup of gauge transforma- 
tions not depending on x"). This has some technical inconvenience but now any periodic 
field can be described in the gauge ( |2.47| ). Furthermore, if we restrict the class of possible 
periodic fields by the condition A_ = 0, disregarding the describe consideration, we come 
to canonical theory with even more complicate constraints if zero modes are taken into 



account 17, 18, U9L 



From the point of view of LF canonical formalism the variables A 1 ! are "coordinates". 
Therefore one can restrict their values by the condition (|2.47| ) directly in the Lagrangian 
without loosing any equations of motion. Let us introduce the denotations 



D_A% = (<9_ - ig(v l - v j ))A%, 



(2.49) 



D^ = {d^-ig^-v^))Al 

where obviously 

D_Af = d_Af. 

Also for any function f(x~) = f(x + ,x~,x x ) periodic in x~ we denote 



(2.50) 



(2.51) 



/(o) = ^ J dx f(x~ 

-L 



(2.52) 



[/(*-)] =/(*-)-/ ( o). 



Obviously, 



(2.53) 



dx~ [f{x-)\ = 0. 



(2.54) 



After gauge fixing ( |2.47| ) the Lagrangian ( 2.45 ) takes the form 



A] 



L= I d 2 x I dx~hj2(d 

-L { i 

+ E(^-4)^ + 2 E[^ 



A\. 



+ 2 E (^-4?) d+A% + E [0- 

i,3, i¥=3 i 



A'_ 



d k d- 



A\ 



+ 



+ 
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+2 £ A^I^-^^ 

+ 1 d 2 J2L]T(d + ^ 2 -4L£(^ 

+i<? / dx- £ (4^-4" - (D_Af ) ) ) . (2.55) 

Here we have ignored some unessential surface terms. 

Variation of the Lagrangian in Af at i ^ j leads to constraints, the solution of 
which can be written in the form 



[A**] = dZ' 



dkd-[Afi-ig E (AfD.A{ t -{D.Af)A{' 



(2.56) 



A 



'j 



(2.57) 



The operator <9_ x is completely defined, as before, by the condition ( |2.28| ) being well 
defined on functions [/(#)]. The operator Dz 1 after Fourier transformation in x~ is 
reduced to the multiplication by (i(k- — g(v l — v 3 ))) . Therefore it has, in general, no 
singularities for any = n(jr/L) with integer n. 

Substituting the expressions ( P-56j ), ( |2.57| ) into the Lagrangian fl2.55|) , we exclude from 
it the quantities [A™] and A % 1 at % ^ j. The variation of the Lagrangian in ^4+(o) leads to 
the constraint 



Q ii (x+,x ± ) = -2 



2Ld k d k v* + tg j dx~ E [A\ 3 D_A{ 1 - (D_A\ 3 )A{'' 
-t, i',i'¥*i 



0. (2.58) 



It is a first class constraint which can be posed on physical state vectors after quantization. 
Therefore we can keep the term with this constraint in the Lagrangian. 

Now we must put in the standard canonical form the terms of the Lagrangian 



d 2 x I dx-{2 y E(d-[A%])d + [A%]+2 £ (D-Ay)d + A% \ . 



(2.59) 



It can be reached by going to Fourier transform 

[At(x-)] = — = £ k - 1/2 {«£(*-) eM-ik-X-) + a\ + {k.) exp^x")} , (2.60) 
2V2L k _ =7T / L 
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where we sum over k- = nn/L, n = 1, 2, . . ., and 



#J 2V2L 



1/2 



f — g(-u 4 — J (&_) exp(zA;_x ) + 



^k—>g(v 1 — v?) 



-1/2 



+ ^2 {k- — g(v J ' — v r ) ) ' a 3 ^(k-) exp(— ik-X ) 

A;_ >g{vi —v l ) 



(2.61) 



where we sum over all k_ = mr/L satisfying corresponding inequalities. The expression 
Q2.59 ) takes the form 



/[ oo 

{ i fc_=7r/L 



al + (A^)5 + oi(*-)) + 



+ E E (a?(*_)9 + ay + (fe_)-ojf + (k_)d + a?(k_ 

i,j,i¥=j k->g(v l -vi) 



Further, in the Lagrangian ( |2.55|) there is a part 



The "momentum" conjugated to v l is 

5L V 



k(0) > 



Hence, 



The corresponding part of the Hamiltonian equals to 

/*■„ = / d 2 x x (V'd+v*) -L V = J d^LY + ^4\o)) 

and the corresponding part of canonical Lagrangian is 

L v = J d 2 x ± ^ (V i d + v i ^j — H v . 



(2.62) 



(2.63) 



(2.64) 



(2.65) 



(2.66) 



(2.67) 



Excluding from the Lagrangian Q2.55| ) the quantities [A™] and A % + (at j / j) vie the 
eqns. (|2.60|) , (|2.61|) , replacing the terms ( |2.59|) by the expression (|2.62|) and the part ( |2.63|) 
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by the expression ( J2.67Q , we obtain the result 

/| oo 
d 2 x L E (aUk^)d + al + (k. 

{ i k-=iv/L 



) - 4 + (k„)d + ai(k-)) + 



+ E E [a{(k-)d + ar (A:_) - ar(k-)d + a%(k-)) + 

+ E^+^+E^+(o)^4-^ 



(2.68) 



where Q M are defined by ( |2.58j ), and the Hamiltonian H = P + is equal to 



H = d 2 x dx- \J2 [9- 



A' 



+ E (n .v!) i) a 1 : - - £ f^f™ \ + 



+2Ljd 2 x ± J2(jl + dkAt {0) )' 



(2.69) 



It is implied that instead of the quantities and A l + (at i ^ j) one uses the 

expressions ( |2.56| ), fl2.57| ) and the A % 1 are expressed in terms of a k , a l k + , 0%, 0% , (at 
i 7^ j) and of A"/ \ with the help of eqns. (|2.60|) , (|2.61|) and 



A\ 



A[ 



+ A 



fc(0)" 



(2.70) 



It is seen from the formulae ( |2.68j ) that a\ + , a^,, , play the role of creation 
and annihilation operators. After quantization they satisfy the following commutation 
relations (at x + = const): 



a\(k-,x x ),ail + {kf_,x , - L ) = ^S kl S k _, k/ _S 2 (x ± -x' ± ) 



(2.71) 



o-k(k-,x ),cli J (k_,x ) 



S ll S"5 kl 5 k _^_5 2 (x ± -x /± ), i^j, «7/. (2.72) 



Also we have 



VHx^),vHx 



-iS ij S 2 (x J 



— x 



(2.73) 



Remaining commutators are equal to zero. 

The operator of the momentum P_, defined by 



P- = d x dx T. 



(2.74) 
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acts on physical states \^f), satisfying the condition 

Q ii (x ± )\^f) = 0, (2.75) 
equivalent to the canonical operator 

P c . an =Jd 2 x(j: E k.al + (k.)al(k.) + £ £ fc_af , (2-76) 

\ i fc_=7r/L k->g(v i —vi) j 

where the normal ordering was made. 

Physical vacuum \ satisfies the relations 

al(k.,x ± )\Q) = 0, (2.77) 



(k_,x ± )\Q) = 0, i^j, 



(2.78) 



and the condition (|2.75| ). 

This scheme is connected with the following essential difficulty. The zero modes 
x L ) are present in the Lagrangian ( [2.68|) and in the Hamiltonian (|2.69| ) but the 

Therefore new constraints arise 



a a 

^(o) 

derivatives d+A 1 ^ are absent there 



5H 



(2.79) 



These constraints are of the 2nd class and they must be solved with respect to A % ^ and 
then the A % ^ have to be excluded from the Hamiltonian. The constraints ( |2.79| ) are very 
complicated and explicit resolution of them is practically impossible. The application of 
Dirac brackets does not simplify this. 

Due to this difficulty a practical calculation usually ignores all zero modes from the 
begining. It makes the approximation worse. It is interesting that in the framework of 
lattice regularization it is possible to overcome the difficulties caused by the constraints 
(|2.79|) This question will be considered in Sect. 5. 



3 Limiting transition from the theory in Lorentz coordinates 
to the theory on the Light Front 

To clarify the connection between the theory in Lorentz coordinates in Hamiltonian form 
and analogous theory on the LF we perform the limiting transition from one to the 
other. Here this transition is considered in the fixed frame of Lorentz coordinates by 
introducing states that move at a speed close to the speed of light in the direction of 
the x 3 axis. Constructing the matrix elements of the Hamiltonian between such states 
and studying the limiting transition to the speed of light (an infinite momentum), we can 
derive information about the Hamiltonian in the light-like coordinates. This information 
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also takes into account the contribution from intermediate states with finite momenta. 
Here, we illustrate the results of such an investigation using (1 + l)-dimensional theory of 
scalar field with the A<^ 4 -interaction. Instead of x 3 we denote analogous space coordinate 
by x 1 . The generalization of the method to (3 + 1) -dimensional Yukawa model is discussed 
briefly at the end of this section. The limiting transition studied here is accomplished 
approximately by subjecting the momenta p\ to an auxiliary cutoff that separates fast 
modes of the fields (with high p x values) from slow modes (with finite p\ values). This 
cutoff is parametrized in terms of the quantities A, A 1; and 5 and the limiting-transition 
parameter 7] (rj > 0, rj — > 0): we have ?7 _1 Ai > \pi\ > r]~ 1 5 for the fast modes and pi < A 
for the slow modes (A ^> 5). For 77 — ^ 0, the inequality i]~ 1 5 > A holds, so that the above 
momentum intervals are separated. The field modes with the momenta r]~ 1 8 > \p\\ > A 
are discarded. This procedure is justified by the fact that the resulting Hamiltonian in 
the limit 77 — > reproduces the canonical light-front Hamiltonian (without zero modes) 
when only the fast modes are taken into account and is consistent with conventional 
Feynman perturbation theory for 5^0. Therefore, even an approximate inclusion of 
the other (slow) modes may provide a description of nonperturbative effects, such as 
vacuum condensates. The effective light-front Hamiltonian obtained here for the model 
under consideration differs from the canonical Hamiltonian only by the presence of the 
vacuum expectation value of the scalar field and by an additional renormalization of the 
mass of this field. The renormalized mass involves the vacuum expectation value of the 
squared slow part of the field. Masses of bound states can be found by solving Schrodinger 
equation 



2 

m 



= — (3.1) 

with obtained Hamiltonian P + . 

We start from the standard expression for the Hamiltonian of scalar field (p(x) in 
(1 + l)-dimensional space-time in Lorentz coordinates x^ = (x ,^ 1 ), at x° = 0: 

H=:J Sx Qn 2 + \{d^Y + ^y. 2 + A^) :, (3.2) 

where II (a; 1 ) are the variables that are canonically conjugate to ^(x 1 ) = ip(x° = 0, x 1 ), and 
the symbol : : of the normal ordering refers to the creation and annihilation operators 
a and a + that diagonalize the free part of the Hamiltonian in the Fock space over the 
corresponding vacuum |0). These operators are given by 

^{x 1 ) = ^r=\ d Pl (m 2 + p 2 )- 1 / 4 [a( Pl ) expHpxrr 1 ) + h.c] , (3.3) 

V 4:71 J 1 J 



n(x 1 ) = / dpi{m 2 + p\) 1 ^ 4 a(pi) exp(— ipix 1 ) — h.c 

V 47T J 1 



(3.4) 



where a(pi)\0) = 0. 
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To investigate the limiting transition to the light-front Hamiltonian (defined at x + = 
0), it is more convenient to go over from the Hamiltonian (|3.2|) to the operator H + P\ = 
\^2P + , where the momentum Pi has the form 

Pi = J dpia + (pi)a(pi)pi. (3.5) 

Applying the above parametrization of high momenta in terms of 77,77 — > 0, to p\ we 
can then consider the transition to an infinitely high momentum of states as a limit of the 
corresponding Lorentz transformation with parameter 77. To be more specific, we have 
Pi — > (—r]^/2)^ 1 q_, where g_ is a finite momentum in the light-like coordinates, and 

lim^^-^^KF + PO^obi)) = (q'_\(P + ) x+=0 \q-). (3.6) 

It follows that the eigenvalues E + of the operator (P + ) X + =Q that correspond to the mo- 
mentum q_ are obtained as the corresponding limit of the eigenvalues E{rj) of the operator 
(H + Pi) x o =0 at momentum px- 

E + = lim(77v / 2) _1 P(77). (3.7) 

In the following, we consider this limiting transition as a part of the eigenvalue problem 
for the operator H + P%, using perturbation theory in the parameter rj. Separating the 
Fourier modes of the field into fast and slow ones, as is indicated above, and neglecting the 
region of intermediate momenta (r]~ 1 5 < \pi\ < r]~ 1 Ai is the region of the fast modes, and 
\pi\ < A is the region of the slow modes), we can substantially simplify this perturbation 
theory. The 77 dependence of the field operators and Hamiltonian can then be determined 
by making, in the region of fast momenta, the change of the variables as 

Pl = 7] -^k, a( Pl ) = ^/rja(k), 5 < \k\ < A l5 [a(k), a + (k')} = 6(k - k'). (3.8) 

The fast part [^(x 1 )]/ of the operator fi^x 1 ) is estimated as 

[<p(x% = 0(y) + 0(t? 2 ), y = r i - 1 x\ (3.9) 

(p{ y ) = (4tt)- 1/2 J dk\k\- 1/2 [a{k) exp(-iky) + h.c.}. (3.10) 

We denote the slow part of the field <p by <p (<p = [(p]f + <p). Substituting formulas 
(|3.8|)-(|3TT(]|) into Hamiltonian (|3.2p, we obtain 



H + P 1 = tT 1 {h + rjhi + V 2 h 2 + •••), (3.11) 



Ai 



h = 2j dkd + {k)d{k)k, (3.12) 
6 
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h 1 = {H + P 1 



'(p=cp,U=U 



(H + Pi) 



(3.13) 



dk 



<5<|fc|<Ai 



m 
2\k\ 



~a + (k)~a(k) + : A dy (p\y) + A<p(O)0%y) + 6<p 2 (0)<p 2 (y) : . (3.14) 



Prior to performing integration with respect to y, we formally expanded the operators 
(fix 1 ) = <p{i]y) in Taylor series in the variable r]y and estimated their orders in the 
parameter t] at fixed y. Such estimates can be justified at least in Feynman perturbation 
theory. The operator (H+Pi in Q3.13j ) is defined in such a way that its minimum eigenvalue 
is zero. Let us consider perturbation theory in the parameter 77 for the equation 



{H + P 1 )\f) = E\f) 



(3.15) 



under the condition that the states |/) have, as in formula (|3.6| ), a negative value of 
Pi proportional to rj" 1 for i] — > and also describe the states with a finite mass. The 
expansions of the quantity E and the vector |/) in power series in 7] can then be written 



as 



E = r,- l Y.V n Em \f) = Y,V n \fn)- 

n=2 n=0 



We arrive at the system of equations 

fco|/o) = 0, 



(3.16) 



(3.17) 



ho\fi) + hi\f ) = 0, 



(3.18) 



h \f2) + h 1 \f 1 ) + (h 2 -E 2 )\f ) = 0, 



(3.19) 



To describe solutions to these equations, we use the basis generated by the fast-field 
operators d + {k) over the vacuum |0) and the slow-field operators ip and n over the vacuum 
\v) that corresponds to the Hamiltonian (H + Pi). The vectors of this basis can be 
symbolically represented as 



a+ ...a + \0)ip...ipf[.. .IL\v). 



(3.20) 



By virtue of ( |3.12j) , the manifold of solutions |/o) to equation ( p. 17Q is reduced to the set 
of vectors ( 3.20| ), which do not contain the operators d + (k) with k > 5. Let Vo be the 
projection operator onto this set. According to equation ( |3.18j ), we then have 



V h \fi) = (H + Pi)\f ) = 0. 



(3.21) 



Equation ( |3.21| ) requires that the vectors |/o) be the lowest eigenstates of the operator 
(H+Pi), that is, linear combinations of basis vectors (|3.20|) including neither the operators 
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a + (k) with k > S nor the operators <p and it. We denote the projection operator on this 
set of vectors ( |3.20|) by V' . To determine the quantity E 2 which we are interested in, it 



is sufficient to consider the Po-projection of equation (|3.19| ). Taking into account (|3.17| ) 



( |3.20| ), and ( |3.21 ), we find that E 2 appears as a solution to the eigenvalue problem 

V' h 2 \f ) = E 2 \f ). (3.22) 



Thus, in accordance with ( |3.6| ) and (p.7|), the operator Vq}i 2 Vq plays the role of the effec- 



tive light-front Hamiltonian P+. Substituting formula ( |3.14|) for the operator h 2 into the 



expression for , we take into account that, between the projection operators V' Q , the 
contribution of the field modes with positive momenta (k > 5) vanishes and that the prod- 
ucts of the operators of the slow part of the field can be replaced with their expectation 
values for the vacuum \v). In addition, we note that, under the Lorentz transformation 
corresponding to the limiting transition 77 — > in formula (|3.6| ), the variable y goes over 
into the light-like coordinate y~ = —y/y/2, the momenta k go over into the light-like 
momenta g_ = —y2k, and the corresponding coordinate y + vanishes at x° = (for finite 
values of y~). Going over to the operators A(q ) = 2~ 1 / 4 a(— k) for k < —5 (q- > 5\/2) 
and to the corresponding field $(y + = 0, y~) = (f-(y), where y5_ is the part of the (p 
containing only the modes with negative momenta (k < —S), we obtain the effective 
Hamiltonian Pi in the form 



jdy- [m 2 + 12A(: 2 :)„] $ 2 + 4\(^) v ^ + A$ 4 | :, (3.23) 



where (. . .)„ is the expectation value for the vacuum \v). Expression Q3.23j ) coincides with 
the canonical effective Hamiltonian if, in the latter, we take into account the shift of the 
field by the constant ((f) v and the change in the mass squared by 12A [(: ip 2 :) v — ((f) 2 ]- 

Analogous results were obtained for Yukawa model in (3 + l)-dimensional space-time 
H . In regularization of Pauli-Villars type, introducing a number of nonphysical fields with 
very large masses. The absence of essential difference between the Hamiltonian obtained 
via limiting transition and canonical LF Hamiltonian is connected with this choice of 
regularization. Other regularizations can lead to more complicated results. 

This method of limiting transition can not be directly expanded to gauge theories, 
because the approximations used for nongauge theories are not justified. 



4 Comparison of Light Front perturbation theory 
with the theory in Lorentz coordinates 

As is already known, canonical quantization in LF, i.e., on the x + = const hypersurface, 
can result in a theory not quite equivalent to the Lorentz-invariant theory (i.e., to the 
standard Feynman formalism). This is due, first of all, to strong singularities at zero 
values of the "light-like" momentum variables Q_ = 4|(Qo _ Qz)- To restore the equiva- 
lence with a Lorentz-covariant theory, one has to add unusual counter-terms to the formal 
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canonical Hamiltonian for the LF, H = P + . These counter-terms can be found by compar- 
ing the perturbation theory based on the canonical LF formalism with Lorentz-covariant 
perturbation theory. This is done in the present section. The light-front Hamiltonian 
thus obtained can then be used in nonperturbative calculations. It is possible, however, 
that perturbation theory does not provide all of the necessary additions to the canonical 
Hamiltonian, as some of these additions can be nonperturbative. In spite of this, it seems 
necessary to examine this problem within the framework of perturbation theory first. 

For practical purposes a stationary noncovariant light-front perturbation theory, which 
is similar to the one applied in nonrelativistic quantum mechanics, is widely used. It 



was found |20| , |21j |22f that the "light-front" Dyson formalism allows this theory to be 
transformed into an equivalent light-front Feynman theory (under an appropriate regu- 
larization). Then, by re-summing the integrands of the Feynman integrals, one can recast 
their form so that they become the same as in the Lorentz-covariant theory. (This is not 
the case for diagrams without external lines, which we do not consider here.) Then, the 
difference between the light-front and Lorentz-covariant approaches that persists is only 
due to the different regularizations and different methods of calculating the Feynman in- 
tegrals (which is important because of the possible absence of their absolute convergence 
in pseudo- Euclidean space). In the present section, we concentrate on the analysis of this 
difference. 

A light-front theory needs not only the standard UV regularization, but also a spe- 
cial regularization of the singularities Q_ =0. In our approach, this regularization (by 
method ( |1 . 1| ) ) eliminates the creation operators a + (Q) and annihilation operators a(Q) 
with \Q l _\ < e from the Fourier expansion of the field operators in the field representa- 
tion. As a result, the integration w.r.t. the corresponding momentum Q_ over the range 
(—00, —e) U (e, 00) is associated with each line before removing the (^-functions. Differ- 
ent propagators are regularized independently, which allows the described re- arrangement 
of the perturbation theory series. On the other hand, this regularization is convenient 
for further nonperturbative numerical calculations with the light-front Hamiltonian, to 
which the necessary counter-terms are added (the "effective" Hamiltonian). We require 
that this Hamiltonian generate a theory equivalent to the Lorentz-covariant theory when 
the regularization is removed. Note that Lorentz- invariant methods of regularization (e.g., 
Pauli-Villars regularization) are far less convenient for numerical calculations and we shall 
only briefly mention them. 

The specific properties of the light-front Feynman formalism manifest themselves only 
in the integration over the variables Q± = -^(Qo ± Q3), while integration over the trans- 
verse momenta Q± = {Qi, Q2} is the same in the light-front and the Lorentz coordinates 
(though it might be nontrivial because it requires regularization and renormalization). 
Therefore, we concentrate on a comparison of diagrams for fixed transverse momenta 
(which is equivalent to a two-dimensional problem). 

In this section, we propose a method that allows one to find the difference (in the limit 
e —>■ 0) between any light-front Feynman integral and the corresponding Lorentz-covariant 
integral without having to calculate them completely. Based on this method, a procedure 
is elaborated for constructing an effective Hamiltonian in LF in any order of perturbation 
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theory. The procedure can be applied to all nongauge field theories, as well as to Abelian 
and non- Abelian gauge theories in the gauge A_ = with the vector meson propagator 



chosen according to the Mandelstam-Leibbrandt prescription p3|, The question of 
whether the additional components of the Hamiltonian that arise can be combined into a 
finite number of counter-terms must be dealt with separately in each particular case. 

Application of this formalism to the Yukawa model makes it possible to obtain the 
effective light-front Hamiltonian in a closed form. The result agrees with the conclusions 



of the work |2~0|, where a comparison was made of the light-front and Lorentz-covariant 



methods via calculating self-energy diagrams in all orders of perturbation theory and other 
diagrams in lowest orders. Conversely, for gauge theories (both Abelian and non- Abelian) , 
it was found that counter-terms of arbitrarily high order in field operators must be added 
to the effective Hamiltonian. This result may turn out to be wrong if the contributions 
to the counter-terms are mutually canceled. This calls for further investigation, but such 
possibility appears to be very unlikely. 

What we have said above does not depreciate the light-front formalism as applied to 
gauge theories. This is because the only requirement concerning the light-front Hamil- 
tonian is that it correctly reproduces all gauge-invariant quantities rather than the off- 
mass-shell Feynman integrals in a given gauge. However, renormalization of the light-front 
Hamiltonian turns out to be a difficult problem and it requires new approaches. We do 
not examine the possibilities of changing the light-front Hamiltonian by introducing new 



nonphysical fields by a method different from the Pauli-Villars regularization |25[ or the 
possibilities of using gauges more general than A_ = with the Mandelstam-Leibbrandt 
propagator. These points also need to be investigated further. 

4.1 Reduction of light-front and Lorentz-covariant Feynman 
integrals to a form convenient for comparison 

Let us examine an arbitrary IPI Feynman diagram. We fix all external momenta and all 
transverse momenta of integration, and integrate only over Q + and Q_: 

F-\\m! u ' d2Qt fiQ%pk) an 
Yli(2Q i +Q i --Mi + i&)' { ' ' 

We assume that all vertices are polynomial and that the propagator has the form 
z(Q) z(Q) Q + 



or (4 2) 

Q 2 -m 2 + i^ {Q 2 - m 2 + iae)(2Q + Q_ + zae) ' ( ' J 

where z{Q) is a polynomial. A propagator of the second type in ( |4.2j ) arises in gauge 



theories in the gauge A_ = if the Mandelstam-Leibbrandt formalism |23|, |J] with the 
vector boson propagator 



Q 2 + iae V" 1 " 2Q + Q_ + ias 



9w 
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is used. In eq. ( J4.1|) either Mf = mf + Q 1 ^ ^ 0, where rrii is the particle mass, or Mf = 0. 

The function / involves the numerators of all propagators and all vertices with the 
necessary (^-functions, that include the external momenta p k (the same expression without 
the ^-functions is a polynomial, which we denote by /). We assume for the diagram F 
and for all of its subdiagrams that the conditions 

w„ < 0, uj + < 0, (4.3) 

hold, where u + is the index of divergence w.r.t. Q + at Q l _ ^ Vz, and iv n is the index 
of divergence in Q + and Q- (simultaneously); Q± = ^(Qo ± Qz)- The diagrams that 
do not meet these conditions should be examined separately for each particular theory 
(their number is usually finite). We seek the difference between the value of integral 
(|4.1|) obtained by the Lorentz-covariant calculation and its value calculated in light-front 
coordinates (light-front calculation). 

In the light-front calculation, one introduces and then removes the light-front cutoff 
\Q-\ > £ > 0: 



f(Q\p k 



Ui^Q'+Qt-Mi + iee) 



where V £ = Hi (( — oo, —e) U (s, oo)). Here (and in the diagram configurations to be defined 
below) we take the limit w.r.t. e, but, generally speaking, this limit may not exist. In 
this case, we assume that we do not take the limit, but take the sum of all nonpositive 
power terms of the Laurent series in e at the zero point. If conditions ( j4.3| ) are satisfied, 
Statement 2 from Appendix I can be used. This results in the equality 

F lf = limlim [T]d& I TW , tJJ^KL r- (4-4) 

From here on, the momenta of the lines Q % are assumed to be expressed in terms of the 
loop momenta q k , B L is a sphere of a radius L in the g^i-space, and L depends on the 
external momenta. Now, using Statement 2 from Appendix I, we obtain 

r r f(O i r) s ) p-^Y. i Q i +-PY, i Q- 2 
Ftf = lim lim lim lim / TT dq\ \ TT dq k JW \ P ± -= — . (4.5) 

Ve 

To reduce the covariant Feynman integral to a form similar to ( |4.4|) , we introduce a 
quantity F: 

F = lim lim lim / Y\d 2 q k ' . = ^. (4.6) 

Let us prove that this quantity coincides with the result of the Lorentz-covariant calcula- 
tion .F cov . To this end, we introduce the a-representation in the Minkowski space of the 
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propagator 

z(Q l 



JM [ e ™iCQ\Qt-M?+i*)+i(Q* + y++Qly7) dai _ 



—iz —i 



Then we substitute (f4.7|) into ( f4.6|) . Due to the exponentials that cut off g+, q k _ and a % 
the integral over these variables is absolutely convergent. Therefore, one can interchange 
the integrations over q+, q k i and a 1 . As a result, we obtain the equality 



where 



F = lim lim lim / TTdctj (p{a i} p s , 7, (3) e _ae S;' 

ae^O /3^0 7^0 J x x 



<^(ai,p s ,7,/3) = (-i) n f x 
f Y[d 2 q k e Y, t [ ia ^ 2 QXQ-- M f)+ i (QXyt+Q-yr)-^Q+ 2 -PQ~ 



(4- 



Vi=0 



(4.9) 



For the Lorentz-covariant calculation in the a- representation satisfying conditions (|4.3| ), 
there is a known expression [2S ] 



F cov = lim / TT v?, 

ae— >0 J 



where 



<9 



(4.10) 



^cov(«i,p s ) = {~i) n f ( J x 



X 



lim /TTdVeS.N^' 

7,0-0 ./ V 



-Mf)+i(Q* + y++Q*_y7)-~,Ql 2 -f3Ql 2 ] 



Vi=0 



(4.11) 



In Appendix 2, it is shown that in (|4.8| ) the limits in 7 and /3 can be interchanged, in 
turn, with the integration over {c^}, and then with / (— i-^j- After that, a comparison 

of relations Q) and ( CT) , ( pT|) , clearly shows that F = F cov . Considering Q 

and using Statement 1 from Appendix 1, we obtain the equality 



Si 



_QL - M? + iae) ' 

Expression ( |4.12| ) differs from ( |4.4| ) only by the range of the integration over q^. 



(4.12) 
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4.2 Reduction of the difference between the light-front 
and Lorentz-covariant Feynman integrals 
to a sum of configurations 

Let us introduce a partition for each line, 




We call a line with integration w.r.t. the momentum Q l _ in the range (— e, e) (before 
removing the 5-functions) a type-1 line, a line with integration in the range (—00, — e) U 
(e, 00) a type-2 line, and a line with integration over the whole range (— 00, 00) a full line. 
In the diagrams, they are denoted as shown in Figs, la, b, and c, respectively. 

9 H- X i 

abode 



Fig. 1: Notation for different types of lines in the diagrams: "a" is a type-1 line, "b" is a 
type-2 line, "c" is a full line, "d" is an e-line, and "e" is a H-line. 



Let us substitute partition ( |4.13| ) into expression ( |4.4j ) for and open the brackets. 



Among the resulting terms, there is F cov (expression (|4.12 )). We call the remaining 



terms "diagram configurations" and denote them by Fj. Then we arrive at the relation 

Fx - F cov =Y,Fj, where 
j 

Fj = lim lim / UdqH f IT dq k , tJ^^X, ( 4 - 14 ) 

and Vj is the region corresponding to the arrangement of full lines and type-1 lines in the 
given configuration. 

Note that before taking the limit in e, Eqs. ( |4.12| ) and ( (4.14Q can be used successfully: 



first, they are applied to a subdiagram and, then, are substituted into the formula for the 
entire diagram. This is admissible because, after the deformation of the contours described 
in the proof of Statement 1 from Appendix 1, the integral over the loop momenta {q+} of 
the subdiagram converges (after integration over the variables {q 1 !} of this subdiagram) 
absolutely and uniformly with respect to the remaining loop momenta {q^_ }. Therefore, 
one can interchange the integrals over and {q 1 !}- 

Thus, the difference between the light-front and Lorentz-covariant calculations of the 
diagram is given by the sum of all of its configurations. A configuration of a diagram is 
the same diagram, but where each line is labeled as a full or type-1 line, provided that at 
least one type-1 line exists. 
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4.3 Behavior of the configuration as e — > 

We assume that all external momenta p s are fixed for the diagram in question and 

E^-VO, (4.15) 

a' 

where the summation is taken over any subset of external momenta; all of these momenta 
are assumed to be directed inward. 

Let us consider an arbitrary configuration. We apply the term "e-line" to all type-1 
lines and those full lines for which integration over Q_ actually does not expand outside 
the domain (—re, re), where r is a finite number (below, we explain when these lines 
appear). The remaining full lines are called Il-lines. In the diagrams, the e-lines and II- 
lines are denoted as shown in Figs. Id and e, respectively. Note that the diagram can be 
drawn with lines "a" and "c" from Fig. 1 (this defines the configuration unambiguously), 
or with lines "d" and "e" (then the configuration is not uniquely defined). 

If among the lines arriving at the vertex only one is full and the others are type-1 
lines, this full line is an e-line by virtue of the momentum conservation at the vertex. The 
remaining full lines form a sub diagram (probably unconnected). By virtue of conditions 
(|4.15 ), there is a connected part to which all of the external lines are attached. All 



of the external lines of the remaining connected parts are e-lines. Consequently, using 
Statement 1 from Appendix 1, we can see that integration over the internal momenta of 
these connected parts can be carried out in a domain of order e in size, i.e., all of their 
internal lines are e-lines. Thus, an arbitrary configuration can be drawn as in Fig. 2 and 
integral Q4.14|) , with the corresponding integration domain, is associated with it. 



n 



£ 



Fig. 2: Form of an arbitrary configuration: II is the connected subdiagram consisting of 
Il-lines, £ is the subdiagram consisting of e-lines and, probably, containing no vertices. 

Let us investigate the behavior of the configuration as e — > 0. From here on, it is 
convenient to represent the propagator as 



where z(Q) = z(Q) or z(Q) 



(4.16) 



Q 2 -m 2 + iae' 2Q_ + iae/Q+ 

rather than as ( f4.2| ). Then, in (O), Mf = mf + Q* ± 2 ^ and the function / is no longer a 
polynomial. If the numerator of the integrand consists of several terms, we consider each 
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term separately (except when the terms arise from expressing the propagator momentum 
Q l _ in terms of loop and external momenta). 

We denote the loop momenta of subdiagram II in Fig. 2 by q l and the others by k m . 
We make following change of integration variables in ( f4.14| ): 



e k r 



(4.17) 



Then, the integration over k™ goes within finite limits independent of e. We denote 
the power of e in the common factor by r (it stems from the volume elements and the 
numerators when the transformation ( [4.1 7[ ) is made). The contribution to r from the 
expression 1/(2(5- + i as /Q+) (Eq. ( |4.16|) ), which is related to the e-line, is equal to -1. 
We divide the domain of integration over k™ and q l + into sectors such that the momenta 
of all full lines Q\ have the same sign within one sector. 

In Statement 1 of Appendix 1, it is shown that for each sector, the contours of inte- 
gration over q l _ and k™ can be bent in such a way that absolute convergence in q\, k™, 
q l _ and k™ takes place. Since, in this case, the momenta Q l _ of II-lines are separated from 
zero by an e- independent constant, the corresponding Il-line- related propagators and fac- 
tors from the vertices can be expanded in a series in e. This expansion commutes with 
integration. 

It is also clear that the denominators of the propagators allow the following estimates 
under an infinite increase in 



1 



2Q+Q- ~ M 2 + ias 



< 



c\Q+\ 
l 

c e \Q A 



for II-lines, 
for e-lines, 



(4.18) 
(4.19) 



Here c and c are e-independent constants. Note that for fixed finite Q + , the estimated 
expressions are bounded as e — > 0. After transformation ([4.17 ) and release of the factor 
- (in accordance with what was said about the contribution to r), the e-line-related 
expression from (|4.16|) becomes 



2Q_ + ias/Q A 



2Q_ + ia}/(Q + e) 



< 



2\Q. 



where a Q + -independent quantity was used for the estimate (this quantity is meaningful 
and does not depend on e because the value of Q- is separated from zero by an e- 
independent constant). 

We integrate first over q l + , k™ within one sector and then over q_, k™ (the latter 
integral converges uniformly in e). Let us examine the convergence of the integral over 
q\, k™ with canceled denominators of the e-lines (which is equivalent to estimating the 
expressions (|4.19 ) by a constant). If it converges, then the initial integral is obviously 
independent of e and the contribution from this sector to the configuration is proportional 
to e T . 
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Let us show that if it diverges with a degree of divergence a, the contribution to the 
initial integral is proportional to e T ~ a up to logarithmic corrections. To this end, we 
divide the domain of integration over q\ , k™ into two regions: U\, which lies inside a 
sphere of radius A/e (A is fixed), and U 2 , which lies outside this sphere (recall that in our 
reasoning, we deal with each sector separately). Now we estimate (|4.18| ) (like (^4.19|) ) in 

terms of — —- — r (which is admissible) and change the integration variables as follows: 

ce\Q + \ 

q l + ^- £ q l + , k™^±k™. (4.20) 

After e is factored out of the numerator and the volume element, the integrand becomes 
independent of e. Thus, the integral converges. 

One can choose such A (independent of e) that the contribution from the domain {7 2 
is smaller in absolute value than the contribution from the domain U\. Consequently, 
the whole integral can he estimated via the integral over the finite domain U\. Now we 
make an inverse replacement in fl4.20|) and estimate ( |4.19| ) by a constant (as above). Since 
the size of the integration domain is A/e and the degree of divergence is a, the integral 
behaves as e~ a (up to logarithmic corrections), q.e.d. This reasoning is valid for each 
sector and, thus, for the configuration as a whole. Obviously, 

a = maxav, (4-21) 

where a r is the subdiagram divergence index and the maximum is taken over all subdia- 
grams D r (including unconnected subdiagrams for which a r is the sum of the divergence 
indices of their connected parts). In the case under consideration, a r = uj\ + v r , where 
v T is the number of internal e-lines in the subdiagram D r . The quantities uo± are the UV 
divergence indices of the subdiagram D T w.r.t. Q±. 

Above, we introduced a quantity r, which is equal to the power of e that stems 
from the numerators and volume elements of the entire configuration. We can write r = 
u) r _ —fi r + u r + rj r , where // is the index of the UV divergence in Q_ of a smaller subdiagram 
(probably, a tree subdiagram or a nonconnected one) consisting of Il-lines entering D r . 
The term rf is the power of e in the common factor, which, during transformation ( |4.17| ), 



stems from the volume elements and numerators of the lines that did not enter D r . (It 
is implied that the integration momenta are chosen in the same way as when calculating 
the divergence indices of D r .) Then, up to logarithmic corrections, we have 

Fj ~ e\ = min(r, uj r _ - u r + - // + rf). (4.22) 

Consequently, for e — > 0, the configuration is equal to zero if a > 0. Relation ( [4.22j ) allows 
all essential configurations to be distinguished. 

4.4 Correction procedure and analysis of counter-terms 

We want to build a corrected light-front Hamiltonian if™ 1 with the cutoff \Q l _\ > e, which 
would generate Green's functions that coincide in the limit e — > with covariant Green's 
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functions within the perturbation theory. We begin with a usual canonical Hamiltonian 
in the light-front coordinates if if with the cutoff \Q % _\ > e. We imply that the integrands 
of the Feynman diagrams derived from this light-front Hamiltonian coincide with the 



covariant integrands after some resummation |20], ^1, 22| . However, a difference may 
arise due to the various methods of doing the integration, e.g., due to different auxiliary 
regularizations. As shown in Sec. W?\ this difference (in the limit e — > 0) is equal to 
the sum of all properly arranged configurations of the diagram. One should add such 
correcting counter-terms to if if, which generates additional "counter-term" diagrams, 
that reproduce nonzero (after taking limit w.r.t. e) configurations of all of the diagrams. 
Were we able to do this, we would obtain the desired H^ or . In fact, we can only show how 
to seek the H£° T that generates the Green's functions coinciding with the covariant ones 
everywhere except the null set in the external momentum space (defined by condition 
(|4.15|) ). However, this restriction is not essential because this possible difference does not 
affect the physical results. 

Our correction procedure is similar to the renormalization procedure. We assume that 
the perturbation theory parameter is the number of loops. We carry out the correction 
by steps: first, we find the counterterms to the Hamiltonian that generate all nonzero 
configurations of the diagrams up to the given order and, then, pass to the next order. 
We take into account that this step involves the counter-term diagrams that arose from 
the counter-terms added to the Hamiltonian for lower orders. Thus, at each step, we 
introduce new correcting counter-terms that generate the difference remaining in this 
order. Let us show how to successfully look for the correcting counter-terms. 

We call a configuration nonzero if it does not vanish as e — > 0. We call a nonzero 
configuration "primary" if n is a tree subdiagram in it (see Fig. 2). Note that for this 
configuration, breaking any n-line results in a violation of conditions ( 4.15|) ; then, the 
resulting diagram is not a configuration. We say that the configuration is changed if all of 
the n-lines in the related integral ( 4.14j ) are expanded in series in e (see the reasoning above 
Eq. ( [4.18[ ) in Sec. |4.3| ) and only those terms that do not vanish in the limit e — > after 
the integration are retained. As mentioned above, developing this series and integration 
are interchangeable operations. Thus, in the limit e — » 0, the changed and unchanged 
configurations coincide. Therefore, we always require that the Hamiltonian counter-terms 
generate changed configurations, as this simplifies the form of the counter-terms. Using 
additional terms in the Hamiltonian, one can generate only counter-term diagrams, which 
are equal to zero for external momenta meeting the condition \p s _ | < e, because with the 
cutoff used, the external lines of the diagrams do not carry momenta with \p s _\ < e. We 
bear this in mind in what follows. 

We seek counter-terms by the induction method. It is clear that, in the first order 
in the number of loops, all nonzero configurations are primary. We add the counter- 
terms that generate them to the Hamiltonian. Now, we examine an arbitrary order of 
perturbation theory. We assume that in lower orders, all nonzero configurations that can 
be derived from the counter-terms, accounting for the above comment, have already been 
generated by the Hamiltonian. 

Let us proceed to the order in question. First, we examine nonzero configurations 
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with only one loop momentum k and a number of momenta q (see the notation above 
Eq. ( |4.17|) ). We break the configuration lines one by one without touching the other lines 
(so that the ends of the broken lines become external lines). The line break may result in 
a structure that is not a configuration (if conditions ( |4.15| ) are violated) ; a line break may 
also result in a zero configuration or in a nonzero configuration. If the first case is realized 
for each broken line, then the initial configuration is primary and it must be generated 
by the counter-terms of the Hamiltonian in the order under consideration. If breaking of 
each line results in either the first or the second case, we call the initial configuration real 
and it must be also generated in this order. 

Assume that breaking a line results in the third case. This means that the resulting 
configuration stems from counter-terms in the lower orders. Then, after restoration of 
the broken line (i.e., after the appropriate integration), it turns out that the counter- 
terms of the lower orders have generated the initial configuration (we take into account 
the comment on successive application of Eq. fl4.14|) ; see the end of Sec. |4.2|) with the 
following distinctions: (i) the broken line (and, probably, some others, if a nonsimply 
connected diagram arises after breaking the line) is not a Il-line but a type-2 line, due 
to the conditions \pt\ > e; (ii) if, after restoration of the broken line, the behavior at 
small e becomes worse (i.e., a decreased), then fewer terms than are necessary for the 
initial configuration were considered in the above-mentioned series in e. We expand 
these arising type-2 lines by formula ( |4.13| ) and obtain a term where all of these lines 
are replaced by Il-lines or other terms where some (or all) of these lines have become 
type-1 lines. In the latter case, one of the momenta q becomes the momentum k. We 
call these terms "repeated parts of the configuration" and analyze them together with 
the configurations that have two momenta k. In the former case, we obtain the initial 
configuration up to distinction (ii). We add a counter-term to the Hamiltonian that 
compensates this distinction (the counter-term diagrams generated by it are called the 
compensating diagrams). 

If there is only one line for which the third case is realized, it turns out that, in the 
given order, it is not necessary to generate the initial configuration by the counter-terms, 
except for the compensating addition and the repeated part that is considered at the next 
step. If there are several lines for which the third case is realized, the initial configuration 
is generated in lower orders more than once. For compensation, it should be generated 
(with the corresponding numerical coefficient and the opposite sign) by the Hamiltonian 
counter-terms in the given order. We call this configuration a secondary one. Next, we 
proceed to examine configurations with two momenta k and so on up to configurations 
with all momenta k, which are primary configurations. 

Thus, the configurations to be generated by the Hamiltonian counter-terms can be 
primary (not only the initial primary configurations but also the repeated parts analogous 
to them, called primary-like), real, compensating, and secondary. If the theory does not 
produce either the loop consisting only of lines with Q + in the numerator (accounting 
for contributions from the vertices) or a line with Q + n in the numerator for n > 1, then 
real configurations are absent because a line without Q + in the numerator can always be 
broken without increasing a (see Eq. (|4.22| )). It is not difficult to demonstrate that if each 
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appearing primary, real, and compensating configuration has only two external line, then 
there are no secondary configurations at all. 

The dependence of the primary configuration on external momenta becomes trivial 
if its degree of divergence a is positive, the maximum in formula ( [4.21 ) is reached on 



the diagram itself, and a = 0. Then, only the first term is taken into account in the 
above-mentioned series. Thus, not all of the Il-line-related propagators and vertex factors 
depend on k™ and they can be pulled out of the sign of the integral w.r.t. {k™} in (|4.14j ). 
We then obtain 



Ff im = limlim / T\dk r , 



f'(k m , P S ) 

X 



m2QXQL-Mf + ia S ) 



f"(k m ) 



where V £ is a domain of order e in size. Let us carry out transformations ( |4.17] ) and ( |4.20| ) • 
For the denominator of the Il-line, we obtain 



2(± £ + EP+)(EP-) - M 2 + *ae 2(£ fc+)(EP-) ' 

Here we neglect terms of order e in the denominator because the singularity at fc+ = is 
integrable under the given conditions for a and everything can be calculated in zero order 
in e at o = 0. Thus, the dependence on external momenta can be completely collected 
into an easily obtained common factor. 



4.5 Application to the Yukawa model 

The Yukawa model involves diagrams that do not satisfy condition ( [4.31 ). These are 
displayed in Figs. 3a and b. We have u n = for diagram "a" and uo + = for diagram 
"b". 




Fig. 3: Yukawa model diagrams that do not meet condition Q4.3|) . 



Nevertheless, these diagrams can be easily included in the general scheme of reasoning. 
To this end, one should subtract the divergent part, independent of external momenta, 
in the integrand of the logarithmically divergent (in two-dimensional space, with fixed 
internal transverse momenta) diagram "a". We obtain an expression with < (i.e., 
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which converges in two-dimensional space) and uo + = 0, as in diagram "b". This means 
that the integral over q + converges only in the sense of the principal value (and it is 
this value of the integral that should be taken in the light-front coordinates to ensure 
agreement with the stationary noncovariant perturbation theory). This value can be 
obtained by distinguishing the g + -even part of the integrand. 

Two approaches are possible. One is to introduce an appropriate regularization in 
transverse momenta and to imply integration over them; then, it is convenient to dis- 
tinguish the part that is even in four-dimensional momenta q. The other is to keep all 
transverse momenta fixed; then, the part that is even in longitudinal momenta q n can be 
released. For the Yukawa theory, we use the first approach. For the transverse regular- 
ization, we use a "smearing" of vertices, which is equivalent to dividing each propagator 
by 1 + Q\ /A_|_ 2 . In four-dimensional space, diagram "a" diverges quadrat ically. Under 
introduction and subsequent removal of the transverse regularization, the divergent part, 
which was previously subtracted from this diagram, acquires the form G\ + C%p\- 

After separating the even part of the regularized expression, we fix all of the transverse 
momenta again. Then it turns out that diagrams "a" and "b" in Fig. 3 meet conditions 
( f4.3| ) and one can show that after all of the operations mentioned, the exponent a (see 
(f4.22| )) does not decrease for any of their configurations. Hence, they can be included in 
the general scheme without any additional corrections. 

Let us first analyze the primary configurations (see the definition in Sec. f4.4|) . In the 
numerators, /c_ appears only in the zero or one power and there are no loops where the 
numerators of all of the lines contain Consequently, one always has r > 0, // < 0, 
and if > (see the definitions in Sec. |4.3|) . Analyzing the properties of the expression 
oj r _ — uo r + for the Yukawa model diagrams, we conclude from ( ^4.22| ) that o > always 
holds. The general form of the nonzero primary configurations with a = is depicted in 
Fig. 4. Note that they are all configurations with two external line. 
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Fig. 4: Nonzero configurations in the Yukawa model: p is the external momentum, and 7 + 
or 7" symbols on the line indicate that the corresponding term is taken in the numerator 
of the propagator. In configuration "b", the part that is proportional to 7 + is taken. 

Further, it is clear that there are no nonzero real configurations (see the comment at the 
end of Sec. (4.4|) , and it can be shown by induction that there are no nonzero compensating 
or secondary configurations either (the definitions are given in Sec. [4.4| also). Thus, only 
primary or primary-like configurations can be nonzero and all of them have the form 
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shown in Fig. 4. It can be shown that their degree of divergence a is positive and the 
maximum in formula ( [4.2 1|) is reached for the diagram itself. Thus, the reasoning above 
and below formula ( f4.23|) applies to them. Then, denoting the configurations displayed 
in Figs. 4a-d by D a - D d , we arrive at the equalities D a = ^-C a , D b = ^-C&, D c = C c 
and D d = Cd, where the expressions C a - Cd depend only on the masses and transverse 
momenta, but not on the external longitudinal momenta, and have a finite limit as e — > 0. 

Now we assume that D a - Dd are not single configurations but are the sums of all con- 
figurations of the same form and that integration over the internal transverse momenta has 
already been carried out, (with the above- described regularization). In four-dimensional 
space, the diagrams D a and D b diverge linearly while D c and D d diverge quadratically. 
Therefore, because of the transverse regularization, the coefficients C c and Cd in the limit 
of removing this regularization take the form Ci+C^p^, where C\ and C2 do not depend on 
the external momenta (neither do C a , Cb))- Thus, to generate all nonzero configurations 
by the light-front Hamiltonian, only the expression 

H c = C lV 2 + C 2 pl <p* + C 3 $^il), (4.24) 

V- 

should be added, where (p and i[) are the boson and fermion fields, respectively, and Ci, 
are the constant coefficients. 

Comparing (4.24) with the initial canonical light-front Hamiltonian, one can easily 



see that the found counter-terms are reduced to a renormalization of various terms of 
the Hamiltonian (in particular the boson mass squared and the fermion mass squared 
without changing the fermion mass itself). The explicit Lorentz invariance is absent, 
which compensates the violation of the Lorentz invariance inherent, in the light-front 
formalism. 

Note that in the framework of the second approach, mentioned at the beginning of this 
section, one can obtain the same results. The only difference is that in two-dimensional 
space, the contributions from the configurations displayed in Fig. 3 would additionally 
depend on external transverse momenta. However, this dependence disappears after inte- 
gration over internal transverse momenta with the introduction and subsequent removal 
of an appropriate regularization. 

In the Pauli Villars regularization, it is easy to verify that the expression uj r _ — uj r + — 
fi r +rf from ( [4.22| ) increases. This is because the number of terms in the numerators of the 
propagator increases. Then, the contribution from the e-lines does not change, while the 
n-lines belonging to D r make zero contribution to uj t _ — u r + and 7] r , but —1 contribution 
to ji r . Since r > 0, this regularization makes it possible to meet the condition a > 
for the configurations that were nonzero (one additional boson field and one additional 
fermion field are enough). Then it turns out that the canonical light-front Hamiltonian 
need not be corrected at all. 
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4.6 Application to gauge theories 



Let us consider a gauge theory (e.g., QED or QCD) in the gauge A_ = 0. The boson 
propagator in the Mandelstam-Leibbrandt prescription has the form 

1 / Q„5+Q + + QJ+Q + \ 
Q2 + igd [9^ 2Q + Q„ +i& )' 

All of the above reasoning was organized such that it could be applied to a theory like this 
(with fixed transverse momenta Q± ^ 0). It turns out that there are nonzero configura- 
tions with arbitrarily large numbers of external lines. An example of such a configuration 
is given in Fig. 5. 




Fig. 5: Nonzero configuration with an arbitrarily large number of external lines in a gauge 
theory. The symbols 7 -1 on the lines and the symbols + or _L by the vertices indicate that 
the corresponding terms 7 + or 7 + are taken in the numerators of propagators and in the 
vertex factors. 



Indeed, using formula ([4.22 ), we can see that for the configuration in Fig. 5, r = 



and, thus, o < 0, i.e., this is a nonzero configuration. It is also clear that introduction of 
the Pauli-Villars regularization does not improve the situation because it does not affect 
r. 

Thus, within the framework of the above-described method for correcting the canonical 
light-front Hamiltonian of the gauge theory, an infinite number of counter-terms must be 
added to the Hamiltonian. Note, however, that the formulated conditions for the vanishing 
of the configuration are sufficient, but, generally speaking, not necessary. Because of this 
and because of the possible cancellation of different configurations after integration w.r.t. 
transverse momenta, the number of necessary counter-terms may be smaller. 

5 Transverse lattice regularization of Gauge Theories 
on the Light Front 

Ultraviolet regularization of nonabelian gauge theories via introduction of space-time 



lattice is widely used in nonperturbative considerations [27|. This regularization is in- 



troduced in a gauge invariant way. For canonical formulation on the LF we use trans- 



verse lattice regularization of Bardeen- Pearson type [28| which allows the action to be 
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polynomial in independent variables. As before we formulate the theory on the interval 
— L < x~ < L assuming periodic boundary conditions for all fields in x~ . Further we 
consider again the U(N) theory of pure gauge fields because this example is technically 
more simple than the SU(N) theory. 

The components of gauge field along continuous coordinates be taken 

without modification and related to the sites of the lattice. Transverse components are 
described with complex N x N matrices M k (x), k = 1,2. Each matrix M k (x) is related 
to the link directed from the site x — e k to the site x. The transverse vector e k connects 
two neighbouring sites on the lattice being directed along the positive axis x k {\e k \ = a): 



M k (x) 

• ^* > axis x k 

X X 



The matrix M£(x) is related to the same link but with opposite direction: 



x — e k x 



axis x k 



In this section for the index k the usual rule of summation on repeated indices is not 
used, and where it is necessary the sign of a sum is indicated. The elements of these 
matrices are considered as independent variables. For any closed directed loop on the 
lattice we can construct the trace of the product of matrices M k (x) sitting on the links 
and order from the right to the left along this loop. For example the expression 

Tr [M 2 (x)M l (x - e 2 )M+(x - ei )M+(x)} 

is related to the loop shown in fig. 6. 

It should be noticed that a trace related to closed loop, consisting of one and the same 
link passed in both directions, is not identically unity because the matrices M k are not 
unitary (see, for example, fig. 7). 

The unitary matrices U(x) of gauge transformations act on the M and M + in the 
following way: 

M k (x) - M' k (x) = U(x)M k (x)U + (x - e k ), (5.1) 

M+{x) -> M' k + {x) = U(x - e k )M+{x)U + {x). (5.2) 
To connect these matrices with usual gauge field of continuum theory we can write 
M k (x) = 1 + gaB k (x) + igaA k (x), B+ = B k , A+ = A k . (5.3) 
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M+(x - ei) 




x -e 1 -e 2 



M 2 (x) 



M 1 (x-e 2 ) x-e 2 



axis x 



axis x 



Fig. 6: 



x - e k 




x 



Fig. 7: 

Here Afe(x) coincide with transverse gauge field components in a — > limit; B k (x) are 
extra fields which should be switched off in the limit. 

An analog of field strength can be defined as follows: 

G+_ = F + ^x) = d + A_{x) - d^A + (x) - ig[A + (x), A_(x)}, 

G± >k (x) = — [d ± M k (x) - ig (A±(x)M k (x) - M k (x)A±(x - e k ))\ , 



ga 

G 12 (x) = -— 2 [M 1 (x)M 2 (x - e ± ) - M 2 (x)M 1 (x - e 2 )] . 

ga 

Under gauge transformation they transform as follows; 

G + _(x) -> G' + _(x) = U(x)G + _(x)U + (x), 
G±, k (x) -> G'± k (x) = U(x)G±, k (x)U + (x - e k ), 
G 12 (x) - G" 12 (z) = U(x)G 12 (x)U + (x - ei - e 2 ), 

We choose a simplest form of the action having correct naive continuum limit: 

L 



S = a 2 Y J jdx + J dx- Tr + £ (G+ fc G_ fc + G + _ k G +k ) - G+G 12 



(5.4) 



(5.5) 



+ S m , (5.6) 
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where the additional term S m gives an infinite mass to extra fields B k in a — > limit: 

L 

dx-^Tr \(Mt(x)M,(x) - iY 



m 2 (a) 



Li 

J2 I dx + f dx'^Tr (M£{x)M k {x)-l)' 



— + -m 2 (a) /" (fx 1 - J dx + J dx~ ^ Tr (b£) , m(a) — + oo. 
After fixing the gauge: 

<9_A_ = 0, A l l(x) = 8 ij v j (x ± ,x + ) 
this action can be written in standard canonical form: 

S = « 2 E / dx+ [ dx ~ i E ^.(aOd+u^a 



+ £^(x)Q>' i (a;)-W(a;) 



+ 



where 

D_M l k j (x) = (d- - igv\x) + - e fc )) M l k j (x), 

D_Ml j+ (x) = (cL + ^(x) - y(x - e fe )) M£' + (x), 
D_F; j _(x) = (&_ - igv\x) + i^'(x)) F| J l(:r), 

(x) play the role of Lagrange multiplier, 

Q ji {x) = 2D_F J + i _(x) + 



+ — E E M2' + (*)D-Mi r (x) - Mf^ T (i + e k )D.Ml\x + e k ) 
9 a j> k L 



are gauge constraints and 



h = E(^- + ^ + <4 + <4) 



(5.7) 



(5.8) 



(5.9) 



(5.10) 



(5.11) 



(5.12) 



is the Hamiltonian density. 

These constraints can be resolved explicitly by expressing F+_ in terms of other vari- 
ables with exception of the zero mode components -^+_(o) which can not be found from 
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this constraint equation. The zero mode Q % fa(x , x + ) of the constraint remains unresolved 
and it is imposed as a condition on physical states: 

Qf 0) (x L ,x + )\* phya )=0. (5.13) 

In order to complete the derivation of the action in canonical form and to extract all in- 
dependent canonical variables we make the Fourier transformation in x~ of the transverse 
fields M k (x) as follows: 

M*t{x) = -f= £ {e( Pn + gv\x)-gv^x-e k ))Mi l {(x ± ,x + ) + 



n=— 00 



+© (-Pn ~ 9v%x) + gv j (x - e k )) M l n \ + (x ± , x + )} x 

X 

where 



-1/2 . 

Pn + gv\x) - gv 3 {x - e k ) e iPnX , (5.14) 



©whJ; ' P - = T' neZ - (5 - 15) 



Then the action is 



x 1 - v i 



+i E E E KT<9+m^ + 2l £ A« (0) gf 0) - w(x) j> . (s.ie) 

" - i,j k i 



Here the H is obtained from TC via substitution of the expression 

Ft = (Ft - S ij Ft( 0) ) + <^f-(o) (5-17) 

where the F+_ — 8^F^_f \ are written in terms of M^ k , M^ k + , by solving the constraints 
(|5.11|) and using eq. ( |5.14j ). The F™_, j remain independent. The G l ( 2 are also expressed 



in terms of M% k , M^ k , v l via the eqns. (|5~4"D , ( |5.14|) . 

We have the following set of canonically conjugated pairs of independent variables: 



[v\ n i = 4La 2 F^ (0) }, 

{Mt ^ + }. (5.18) 



In the quantum theory these variables become operators which satisfy the usual canonical 
commutation relations: 

[v l (x),Uj(x')] x + =i5i j 8 x x <x ,±, 
[Mg(x), M% + (x% + = 8^ 8"' 8 nn ,8 kk ,8 x x x ,^ (5.19) 
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In this formulation there are no most complicated constraints like ( |2.79| ) for the zero 
modes of the transversal field components. If one goes to the limit a — > 0, these constraints 
reappear in a form which contains quantum operators in a definite order. This order was 
not clear earlier. 

The operators Q l / Q Jx , x + ) have the following form in terms of canonical variables: 

2LQ™ - ) (x ± , x + ) = 

= -AEEEK^ + + e *) - ^(x)) M^{x + e k )M»(x + e k )- 

/a n j k 

- e (p n + gv\x) - gv\x - e k )) M l n { + (x) M l n {(x)} , (5.20) 

where 

One can easily construct canonical operator of translations in the x~ : 

pcan. = £ £ £ 5> n£ (ft. + gv\x) - gv* (x - e„j) M l n i + (x)M l n i(x). (5.22) 

x 1 - n i,j k 

This expression differs from the physical gauge invariant momentum operator P_ by 
a term proportional to the constraint. The operator P is 



P - = a2 E E / dx ~ Tl { G -k G -k) = P- an ' + 2La 2 E v *Q\ 

I 1„ * I 



'(0) 

k _" L X -L I 



= EEEEK + ^)-^'(^- e& ) M^ + (x)M^(x). (5.23) 

x 1 - n i,j k 

Absolute minimum of P_ corresponds to states \v) defined by the conditions: 

Mg(x)\v) = 0, v\x)\v) = v\x)\v). (5.24) 

The Fock spaces constructed over these states by application of creation operators 
M^ k form the full Fock space of this LF formulation. Operators Q % fa are diagonal in this 
Fock space. In the basis 

\uUUU{Ki(^ k{X) \v)\, (5.25) 

where m^ k [x) are nonnegative integer numbers, we can write the constraint equation 
(|5.13|) in the form: 



EEE [z (Pn + gv 3 (x + e k ) - gv l (x)) m% k (x + e k )- 



j k 



e 



(p n + gv\x) - gv ] (x - e k )J m% k (x) =0. (5.26) 
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One can find the eigenvalue p- of the momentum P_ for such basis state: 



P- = 2 £253 \Vn + gv\x) -gv 3 (x-e k ) m% k (x) (5.27) 



'■j 



and require that this value be finite. For the physical states with p_ = we can find a 
state giving minimum to the Hamiltonian. We consider this state as the vacuum state. 
Here the problem is to solve the reduced Schroedinger equation with variables v % and II* 
which are independent of x_ . 

Analog of this lattice formulation in (2 + l)-dimensions with the Hamiltonian modified 
by addition of some effective interaction terms, while ignoring all zero modes in X-, was 
considered in |3(| in the framework of color dielectric model approach j32|. The 



results of mass spectrum calculation at strong coupling can be fitted to the spectrum 
obtained with usual Wilson-Polyakov lattice theory ||27|| . However the continuum limit 
cannot be reached in this model. 

In our lattice formulation the difficulty with going to weak coupling limit (g — > 0) 
becomes more transparent. The resulting Hamiltonian, written in the canonical variables, 
explicitly contains the terms with zero modes (in x_) having coupling constant in the 
denominator, so that usual perturbative limit does not exist. For example the vacuum 
state corresponding to zero vacuum expectation values of M k (x) does not coincide with 
continuum theory vacuum where these values must be unity. A possible solution of this 
problem can be found in a modification of canonical LF formulation excluding zero modes 
in X- while modifying the LF Hamiltonian so that the theory remains equivalent to original 
Lorentz and gauge invariant formulation after removing the regularizations. One can tried 
to do this by comparison of corresponding perturbation theory series. 



6 Conclusion 

The LF Hamiltonian approach to Quantum Field Theory briefly reviewed here is an at- 
tempt to apply a beautiful idea of Fock space representation for quantum field nonpertur- 
batively in the framework of canonical formulation on the LF. The problem of describing 
the physical vacuum state becomes formally trivial in this approach because such vacuum 
state coincides with mathematical vacuum of LF Fock space. 

However the breakdown of Lorentz and gauge symmetries due to regularizations gen- 
erates difficulties in proving the equivalence of LF formalism and usual one in Lorentz 
coordinates. This problem can be solved for nongauge theories but turns out to be very 
difficult for gauge theories in special (LF) gauge which is needed here. Nevertheless we 
hope that these difficulties can be overcome by finding a modified form of canonical LF 
Hamiltonian which generates the perturbation theory equivalent to usual covariant and 
gauge invariant one. 
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Appendix 1 



Statement 1. If conditions t\4-3j) are satisfied, then, for fixed external momenta p s and 
p s _ 7^ Ws, the equality 



f(Q\p s ) e -^ Q + 2 -^ Q - 2 



/ioWs u>/iz7e £/&e expressions appearing in (\A.1.1\) exist and the integral over {q+} on the 



right-hand side is absolutely convergent. It is assumed that the momenta of lines Q % are 
expressed in terms of loop momenta q k , V £ is the domain corresponding to the presence of 
full lines, type-1 lines, and type-2 lines (the definitions are given following formula j\4-13j )), 
Bl is the sphere of radius L, where L > S max \p s _\, and S is a number depending on the 
diagram structure. 

Let us prove the statement. For each type-1 line in ( A.l.l ), we perform the following 
partitioning: 



J dQ l _ = J dQl + (-1) I J dQt + J dQl 



Then both sides of Eq. ([A.l.lj) become the sum of expressions of the same form in which, 
however, the domain V e corresponds to the presence of only full and type-2 lines. It is 
clear that by proving the statement for this V E : (which is done below), we prove the 
original statement as well. 

Let B be a domain such that the surfaces on which Q l _ = are not tangent to the 
boundary B. First, we prove that in the expression 

v £ nB 

the integral over {<?+} is absolutely convergent (here the integral over {q 1 !} is finite because 
X33 > 0, (3 > 0). This becomes obvious (considering conditions ( f4.3|) and the fact that, in 
type-2 lines, the momentum Q % _ is separated from zero) if the contours of the integration 
over {g^} can be deformed in such a way that the momenta Q % _ of the full lines are 
separated from zero by a finite quantity (within the domain V e C\B). In this case, we can 
repeat the well-known Weinberg reasoning f33| . What can prevent deformation is either 
a "clamping" of the contour or the point Q l _ = falling on the integration boundary. 

Let us investigate the first alternative. We divide the domain of integration over q+ 
into sectors such that the momenta of all full lines Q l + have a constant sign within one 
sector. Let us examine one sector. We take a set of full lines whose Q l _ may simultaneously 
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vanish. In the vicinity of the point where Q l _ from this set vanish simultaneously, we bend 
the contours of the integration over {g^ } such that these contours pass through the points 
Q % _ = iB % and the momenta Q % _ of the type-2 lines do not change. Let B l be such that 
B l Q\ > for the lines from the set (for Q l + from the sector under consideration). It is 
easy to check that this bending is possible. (Since the contours of integration over q_ 
are bent and Q % _ are expressed in terms of g*, one should only check that such b k exist, 
where the necessary B % are expressed in the same way, i.e., that B l obey the conservation 
laws and flow only along the full lines). With this bending, rather small in relation to the 
deviation and the size of the deviation region, the contours do not pass through the poles 
because, for the denominator of each line from the set in question, we have 

(2Q i + Q i _ — M} + *ae) -> (2Q\ (q*1 + iB l ) - Mf + *ae) , Q\B' > 0, 

and for the other denominators, the bending takes place in a region separated from the 
point where the corresponding momenta Q l _ are equal to 0. Repeating the reasoning for 
all sets, we can see that there is no contour "clamping". 

The other alternative is excluded by the above condition for B. To make this clear, one 
should introduce such coordinates £ a in the g fc -space that the boundary of the domain B 
is determined by the equation £ l = a = const and then argue as above for the coordinates 
£ Q with a > 2. 

After bending the contours, integral ( A.1.2 ) is absolutely convergent in q+, q k _ if the 
integration in g+ is carried out within the sector under consideration. On pointing out 



that the result, of internal integration in (|A.1.2 ) does not depend on the bending, we add 



the integrals over all sectors and conclude that ( A.1.2 ) converges in {g+} absolutely. 
Now let us prove that if B is a quite small, finite vicinity of the point {g*:} that 



lies outside the sphere Bl, then expression ( |A.1.2|) is equal to zero. We consider the 



momentum Q l _ of one line. Flowing along the diagram, it can ramify or it can merge 
with other momenta. Clearly, two situations are possible: either it flows away completely 
through external lines, or, probably, after long wandering, part of it, Q_, makes a complete 
loop. The former situation is possible only if \Ql_\ < J2r \P-\i where all external momenta 
leaving the diagram (but not entering it) are summed. Obviously, S can be chosen such 
that for {g^} from B, a line exists whose momentum violates this condition. 

The latter situation results in the existence of a loop, where the inequality Q % _ > Q- 
holds for all momenta of its lines and the positive direction of the momenta is along the 
loop. Then the integral over tn of the loop in question can be interchanged with the 
integrals over {g^} (because it is absolutely and uniformly convergent for all g*) and the 
residue formula can be used to perform this integration. Since, for the loop in question, 
the momenta Q l _ of the lines of this loop are separated from zero and are of the same 
sign, the result is zero. This has a simple physical meaning. If we pass to stationary 
noncovariant perturbation theory, we find that only quanta with positive Q_ can exist. 
In this case, external particles with positive p_ are incoming and those with negative p_ 
are outgoing. Then, the momentum conservation law favors the occurrence of the first 
situation. 
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The entire outside space for B can be composed of the above domains Bl (everything 
converges well at infinity due to the factor exp(— f3J2i Q- ))• Thus, on the left-hand side 
of flA.l.lp , one can substitute the integration domain V £ fl Bl for V e , set the limit in 7 
under the sign of integration over {q%} because of its absolute convergence, and also set 
the limit in f3 under the integration sign because the domain of the integration over {q 1 !} 
is bounded. Thus, we obtain the right-hand side. The statement is proved. 

Statement 2. IfV e corresponds to the presence of type- 2 lines alone, then, under the 
same conditions as in Statement 1, the equality 

KQ\p s ) 

r2 



e 



_Q l _ -Mf+ iae) 



V e nB L 

is valid. 

The proof of this statement is analogous to the second part of the proof of Statement 1. 



Appendix 2 



Statement. If conditions ( \j.3j ) are satisfied, the limits in 7 and (3 in t \4-8j ) can be inter- 
changed (in turn) with the sign of the integral over {ctj} and then with f ( — ^^") ■ 

To prove this, we define the vectors {q+, gi, . . . , q l + , q l _} = S, 
{<3+,Qi, • • • ,Q+,Q1} = [iS + P, and {yt,yi, . . .,y+,y~} = Y, where the vector P is 



built only from external momenta and [i is an I x n matrix of rank I, [i^k-i = Mat 1 = 



t^2k = ^Ik-i- Next, we introduce the following notation: 



A; 



7 \ 

-ieti (3 J ' 



A = diag{Ai, . . . , A n }, A = [i A[i, 



B = [SAP - -i^Y, C = -P l AP + iY l P - i £ c^Mf . 



Then it follows from (|4.9| ) that 



J d 2l S e StAS-2Bt S+ C 



Vi=0 



4 yj (-iJL) e BtA-i B+ c^ = 
! V %/ v / detA 



y l =o 



(A.2.1) 



The function / is a polynomial and we consider each of its terms separately. Up to a 
factor, each term has the form -r^— . . . -J?—. These derivatives act on C and B. The 

' oy il dy ir 

action on C results in the constant factor iN f P, the action on B results in the factor 
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-(l/2)iN t /j,A- 1 B or -(1/4)^/^4" (the latter is the result of the action of two 

derivatives; N, N±, and N 2 are constant vectors). 

It is necessary to prove the correctness of the following three procedures: (i) setting 
the limit in 7 under the integral sign for fixed f3 > 0; (ii) setting the limit in (3 for 7 = 0; 
(iii) setting the limits in 7 and (3 under the signs of differentiation with respect to Y. In 
cases (i) and (ii), one must obtain the bounds 

\p(a u p s j7j (3)\ <p'(a h p s ,P) : (A2.2) 

\p(a llP s ,0,(3)\ <p"(a uP s ), (A2.3) 

where <p' and ip" are functions integrable (for tp f if > 0) in any finite domain over «j, 
with cti > 0. Then, for case (i), we have 

i.e., a limit on the integrated function arises, and, thus, the limit in 7 can be put under 



the integral sign. The situation is similar for case (ii). It is evident from (|A.2.1| ) that the 
function (p can be singular only if the eigenvalues of matrix A become zero. On finding 
the lower bound of these eigenvalues, one can prove through rather long reasoning that 
bounds fljOgfl , (TCT3Q exist if condition Q is satisfied. 



After the limits in 7 and (3 are put under the integral sign, it is not difficult to 
interchange them with the differentiation with respect to Y . One need do it only for 
on > (for each i) and, in this case, one can show that the eigenvalues of the matrix A 
are nonzero and (p is not singular. 
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